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Abstract 

We study in detail the power of the reactor experiment KamLAND for discrim- 
inating existing solutions to the SNP and giving accurate information on neutrino 
masses and mixing angles. Assuming the expected signal corresponding to various 
"benchmark" points in the 2 dimensional (Am 2 , tan 2 9) mixing plane, we develop a 
full-fledged x 2 analysis which include KamLAND spectrum and all the existing solar 
evidence. A complete modelling of statistical and known systematics errors for 1 and 
3 years of data taking is included, exclusion plots are presented. 

We find a much higher sensitivity especially for values of Am 2 lying in the central 
part of the LMA region. The situation would be more complicate for values closer 
to the border of the LMA region (the so called HLMA region, i.e. Am 2 < 2 x 10 5 
and Am 2 > 8 — 9 x 10~ 5 or tan 2 6 far from ~ 0.5). In this case kamLAND, with or 
without solar evidence, will only be able to select multiple regions in the parameter 
space, in the sense that different possible values of the parameters would produce the 
same signal. 
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1 Introduction 



The publication of the recent SNO results |1|-|3| has made an important breakthrough 
towards the solution of the long standing solar neutrino [14, 3^-42| problem (SNP) 
possible. These results provide the strongest evidence so far for flavor oscillation in 
the neutral lepton sector. 

In the near future the reactor experiment KamLAND @,[5| is expected to further 
improve our knowledge of neutrino mixing. In fact it should be able to sound the 
region of the mixing parameter space corresponding to the so called Large Mixing 
Angle (LMA) solution of the solar neutrino problem (Am 2 ~ 10 -5 — 10 -4 eV 2 and 
tan 2 9 ~ 10 -1 — 1) more profoundly. This is of prime interest due to the fact that 
the LMA region in the one preferred by the solar neutrino data at present. 

The previous generation of reactor experiments (CHOOZ ||, Palo Verde |8|), 
performed with a baseline of about 1 km. They have attained a sensitivity of 
Am 2 < 10~ 3 eV 2 |],|7|] and, not finding any dissapearence of the initial flux, they 
demonstrated that the atmospheric neutrino anomaly |ll| is not due to muon-electron 
neutrino oscillations. The KamLAND experiment is the successor of such experi- 
ments at a much larger scale in terms of baseline distance and total incident flux. 

This experiment relies upon a 1 kton liquid scintillator detector located at the old, 
enlarged, Kamiokande site. It searches for the oscillation of antineutrinos emitted 
by several nuclear power plants in Japan. The nearby 16 (of a total of 51) nuclear 
power stations deliver a v e flux of 1.3 x 10 6 cm _2 s _1 for neutrino energies E v > 1.8 
MeV at the detector position. About 78% of this flux comes from 6 reactors forming 
a well defined baseline of 139-214 km. Thus, the flight range is limited in spite of 
using several reactors, because of this fact the sensitivity of KamLAND will increase 
by nearly two orders of magnitude compared to previous reactor experiments. 

It has been estimated that even in the most conservative scenario Q, in which 
the background has to be determined from the reactor power fluctuations, the LMA 
solution is completely within the estimated sensitivity after 3 years of data taking. 
Moreover, KamLAND should be able @ to determine the mixing angle and mass dif- 
ference with a 20% accuracy at 99% of confidence level (CL). As has been underlined 
in [^] and later in [13], the problems might arise if the value of the squared mass 
difference (Am 2 ) lies in the upper region of the LMA solution, the so called HLMA 
region. 

Let us briefly recall some model independent conclusions obtained from the results 
of SNO 111. Different quantities can be defined in order to make the evidence for 
disappearance and appearance of the neutrino flavors explicit. Letting alone the SNO 
data, from the three fluxes measured by SNO is possible to define two useful ratios: 
4>ccl4>ESi 4>ccI4>nc-, deviations of these ratios with respect to their standard value 
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are powerful tests for occurrence of new physics. For the first ratio, one obtains [14] 

^cc _ n 7 o+o.io 
4>bs ~ - 007 ' 

a value which is ~ 2.7 a away from the no-oscillation expectation value of one. The 
ratio of CC and NC fluxes gives the fraction of electron neutrinos remaining in the 
solar neutrino beam, the value obtained in Ref. [14]: 

0.34+f 



- U.d4_ . 04 , 



this value is far away from the standard model case. 

Finally, if in addition to SNO data one consider the flux predicted by the solar 
standard mode one can define, following Ref. [16], the quantity sin 2 a, the fraction of 
" neutrinos which oscillated into active ones", one finds the following result: 

2 _ $NC ~ $CC _ Q 92+0 .39 

sm a — — u.yz_ 20 , 

^SSM — ^CC 

where the SSM flux is taken as the 8 B flux predicted in Ref. |17[ . The central value 
is clearly below one (only-active oscillations). Although electron neutrinos are still 
allowed to oscillate into sterile neutrinos the hypothesis of transitions to only sterile 
neutrinos is rejected at nearly 5a. 

The aim of this work is to study the KamLAND discriminating power, to un- 
derstand in which regions of the parameter space still allowed by the solar neutrino 
experiments KamLAND might give satisfactory accuracy. The structure of this work 
is the following. In section |2| we discuss the main features of KamLAND experi- 
ment that are relevant for our analysis: we derive updated numerical expressions for 
the reactor fuel cycle-averaged antineutrino flux and the absortion antineutrino cross 
section. The next section is devoted to the salient aspects of the procedure we are 
adopting. The results of our analysis are presented and discussed in section || and, 
finally, in section || we draw our conclusions and discuss possible future scenarios. 



2 A Kamland overview 

Electron antineutrinos from nuclear reactors with energies above 1.8 MeV are mea- 
sured in KamLAND by detecting the inverse /3-decay reaction V e + p — > n + e + . The 
time coincidence, the space correlation and the energy balance between the positron 
signal and the 2.2 MeV 7-ray produced by the capture of a already-thermalized neu- 
tron on a free proton make it possible to identify this reaction unambiguously, even 
in the presence of a rather large background. 

The two principal ingredients in the calculation of the expected signal in Kam- 
LAND are the reactor flux and the antineutrino cross section on protons. These 
ingredients are considered below. 
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2.1 The reactor antineutrino flux 



We first describe the flux of antineutrinos coming from the power reactors. A number 



of short baseline experiments ( Ref.[18] and references therein) have measured the 
energy spectrum of reactors at distances where oscillatory effects have been shown 
to be inexistent. They have shown that the theoretical neutrino flux predictions are 
reliable within 2% f|. 

The effective flux of antineutrinos released by the nuclear plants is a rather well 
understood function of the thermal power of the reactor and the amount of thermal 
power emitted during the fission of a given nucleus, which gives the total amount, and 
the isotopic composition of the reactor fuel which gives the spectral shape. Detailed 
tables for these magnitudes can be found in Ref. ||l"8|] . 

For a given isotope (j) the energy spectrum can be parametrized by the following 



expression [19] 



dNl/dE v = exp(a + a x E v + a 2 E 2 v ) (1) 

where the coefficients a% depend on the nature of the fissionable isotope (see Ref. fill 
for explicit values). Along the year, between periods of refueling, the total effective 
flux changes with time as the fuel is expended and the isotope relative composition 
varies. The overall spectrum is at a given time 

dN v \ - dNj 
dE v ~ . 2-* C ^>dE u - 

j=isotopes 

To compute a fuel-cycle averaged spectrum we have made use of the typical evolution 
of the relative abundances c,-, which can be seen in Fig. 2 of Ref.[18|. This averaged 



spectrum can be again fitted very well by the same functional expression (|T]). The 
isotopic energy yield is properly taken into account. As the result of this fit, we 
obtain the following values which are the ones to be used in the rest of this work: 

a = 0.916, oi = -0.202, a 2 = -0.088. 

Although individual variations of the Cj along the fuel cycle can be very high, 
the variation of the two most important ones is highly correlated: the coefficient 
increases in the range ~ 0.5 - 0.7 while c( 239 Pu) decreases ~ 0.4 - 0.2. This 
correlation makes the effective description of the total spectrum by a single expression 
Eq.Q useful. With the fitted coefficients a, above, the difference between this effective 
spectrum and the real one is typically 2 — 4% along the yearly fuel cycle. 

2.2 Antineutrino cross sections 

We now consider the cross sections for antineutrinos on protons. We will sketch 
the form of the well known differential expression and more importantly we will 
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give updated numerical values for the transition matrix elements which appear as 
coefficients . 

In the limit of infinite nucleon mass, the cross section for the reaction V e + p — > n + e~ 
is given by [HJ a(E v ) = kE e +p e + where E,p are the positron energy and momen- 
tum and k a transition matrix element which will be considered below. The positron 
spectrum is monoenergetic: Ep and E e + are related by: E^ = E + + AM, where 
M n , M p are the neutron and proton masses and AM = M n — M p ~ 1.293 MeV. 

Nucleon recoil corrections are potentially important in relating the positron and 
antineutrino energies in order to evaluate the antineutrino flux. Because the an- 
tineutrino flux §{E V ) would typically decrease quite rapidly with energy, the lack of 
adequate corrections will systematically overestimate the positron yield. 

At highest orders, the positron spectrum is not monoenergetic and one has to 
integrate over the positron angular distribution to obtain the positron yield. The 
differential cross section at first order 1/M p is of the form: 



(—) 

Vdcosfl J 



(1) 



<70 

2 L 
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Complete expressions and notation can be found in Ref. [19|. Here we only want to 
pay attention to the overall coefficient gq which is related to the transition matrix 
element k above. 

The matrix transition element can be written in terms of measurable quantities 

as 

k = 2vr 2 log2/ (ml f t 1/2 ) 

where appears the free neutron decay ti/2, the phase-space factor f and m e . The value 
of / = 1.71465 ±0.00015 follows from calculation @, while i 1/2 = 613.9 ±0.55 sec is 
the latest published value for the neutron half-life [23]. This value has a significantly 



smaller error than previously quoted measurements. 

From the values above, we obtain the extremely precise value: 

k = (9.5305 ± 0.0085) x 10" 44 cm 2 /MeV 2 . 

From here the coefficient which appears in the differential cross section is obtained as 
(vector and axial vector couplings / = l,g = 1.26): k = <7o(/ 2 + % 2 )- I n summary, 
the differential cross section which appear in KamLAND are very well known, its 
theoretical errors are negligible if updated values are employed. 



3 The computation of the expected signals 

In order to obtain the expected number of events at KamLAND, we sum the expec- 
tations for all the relevant reactor sources weighting each source by its power and 
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distance to the detector (table II in Ref. |l8| ), assuming the same spectrum orig- 
inated from each reactor. The average number of positrons iVj which are detected 
per visible energy bin AE{ is given by the convolution of different quantities: P, the 
oscillation probability averaged over the distance and power of the different reactors, 
the antineutrino capture cross section given as before, the antineutrino flux spectrum 
given by expression p], the relative reactor-reactor power normalization which is in- 
cluded in the definition of P and the energy resolution of KamLAND which is rather 
good H), we use in our analysis the expression a(E)/E ~ 10%/y/E. 

For one year of running with the 600 ton fiducial mass and for standard nuclear 
plant power and fuel schedule: we assume all the reactors operated at ~ 80% of their 
maximum capacity and an averaged, time-independent, fuel composition equal for 
each detector, the experiment expects about 550 antineutrino events (this number 
agrees with other estimations, i.e. Ref. ||). We will consider this number as our 
KamLAND year in what follows. We will not add any background events as we 
will suppose that they can be distinguished from the signal with sufficiently high 
efficiency. This should be taken with caution, we will dedicate a few words to the 
experimental background at the end of section. 

We compute the expected signal in a set of 0.5 MeV width total positron energy 
bins in the range 2.0 MeV - 8 MeV. Quantitavely, the main information content 
of the shape of the observable spectrum is summarized by the first moment of the 
distribution, the average spectrum energy. This first moment is defined as 

Y^.EiRi (6, Am 2 ) 

— t h 

where Ei is the center positron energy in the i energy bin and the normalized 
signals relative to the non-oscillation case: Ri = Ni{6 , Am 2 ) / 'Nf . 

The expected value of this first moment as a function of Am 2 for some selected 
values of the mixing angle is represented in Fig. [l|(right). We want to illustrate the 
potentiality of this magnitude as an indicator of neutrino oscillations at KamLAND. 
In this plot we graphically see the fact that KamLAND is sensitive to neutrino 
oscillations in the LMA region and in particular to the Am 2 parameter. The first 
moment can vary up ~ 12% in sizeable regions of the parameter space (this can be 
beautifully seen in a 3D plot). A variation of a fraction of this size should be clearly 
identifiable from the KamLAND data after 1 year of effective running. In order to 
see whether this result depends on the choice of bin-size, we have also reproduced it 
with larger, 1 MeV bins, and found that there is no significant difference. 

In Fig.|l] (left) we show the visible positron energy spectrum at KamLAND for 
some chosen oscillation parameters (see Table [j]). We present the integrated signal 
at every 0.5 MeV bin normalized to the non-oscillation expectation. We can see from 
the plot how the shape of the signal is very sensitive to the oscillation parameters. It 
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can greatly change through the LMA region. From the point of view of background 
reduction, in some favorable cases the spectrum is peaked above 5 MeV, this suggests 
the extension of the fiducial energy thresholds beyond the 8 MeV level. 

In addition to the antineutrino signal, two classes of background can be distin- 
guished 

The random coincidence background is due to the contamination of the detector 
scintillator by U, Th and Rn. From MC studies and assuming that an adequate level 
of purification can be obtained, the background coming from this source is expected 
to be ~ 0.15 events/d/kt which is equivalent to a signal to background ratio of ~ 1%. 
Other works [^] conservatively estimate a 5% level for this ratio. More importantly 
for what it follows, one expects that the random coincidence backgrounds will be a 
relatively steeply falling function of energy. The assumption of no background should 
be relatively safe only at high energies (above ~ 5 MeV). 

The second source of background, the so called correlated background is domi- 
nantly caused by cosmic ray muons and neutrons. The KamLAND's depth is the 
main tool to suppress those backgrounds. MC methods estimate a correlated back- 
ground of around 0.05 events/day/kt distributed over all the energy range up to ~ 20 
MeV. 

We will also need the expected signals in the different solar neutrino experiments. 
These are obtained by convoluting solar neutrino fluxes, sun and earth oscillation 
probabilities, neutrino cross sections and detector energy response functions. We 
closely follow the same methods already well explained in previous works [14, 24-26|, 
we will mention here only a few aspects of this computation. We determine the 
neutrino oscillation probabilities using the standard methods found in literature [27], 
as explained in detail in |24| and in [O]. We use a thoroughly numerical method 
to calculate the neutrino evolution equations in the presence of matter for all the 
parameter space. For the solar neutrino case the calculation is split in three steps, 
corresponding to the neutrino propagation inside the Sun, in the vacuum (where the 
propagation is computed analytically) and in the Earth. We average over the neutrino 
production point inside the Sun and we take the electron number density n e in the 
Sun by the BPB2001 model 17 j. The averaging over the annual variation of the orbit 
is also exactly performed. To take the Earth matter effects into account, we adopt 
a spherical model of the Earth density and chemical composition |p8[ . The gluing of 
the neutrino propagation in the three different regions is performed exactly using an 
evolution operator formalism [27]. The final survival probabilities are obtained from 
the corresponding (non-pure) density matrices built from the evolution operators in 
each of these three regions. 

In this analysis in addition to night probabilities we will need the partial night 
probabilities corresponding to the 6 zenith angle bin data presented by SK [43]. 
They are obtained using appropriate weights which depend on the neutrino impact 
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parameter and the sagitta distance from neutrino trajectory to the Earth's center, 
for each detector's geographical location. 



4 Analysis and Results 

In order to study the potentiality of KamLAND for resolving the neutrino oscillation 
parameter space, we have developed two kind of analysis. In the first case (Analysis A 
below) we will deal with the KamLAND expected global signal. We will asumme that 
the experiment measure a certain global signal with given statistical and systematic 
error after some period of data taking (1 or 3 yrs) and will perform a complete x 2 
statistical analysis including in addition the up-date solar evidence. In the second 
case, Analysis B, we will include the full KamLAND spectrum information. Instead 
of giving arbitrary values to the different bins, we will assume a number of oscillation 
models characterized by their mixing parameters (Am 2 , 6). After including the solar 
evidence we will perform the same x 2 analysis as before. 



4.1 Analysis A 

The total x 2 value is given by the sum of two distinct contributions, that is the 
one coming from all the solar neutrino data and the contribution of the KamLAND 
experiment: 



2 _ 2 | 2 
X — Xq T Xglob,KLi 



with 



2 f R e ^-R teo (Am 2 ,e) Y 

X g lob,KL — I ~ I • \ 6 ) 

\ O stat+sys / 

For the "experimental" signal ratio R ex P we assume different values varying from a 
very strong suppression R ~ 0.3 to unobservation of neutrino oscillations R ~ 1. The 
total error a is computed as a sum of assumed systematics deviations, a sys /S ~ 5%, 
mainly coming from flux uncertainty (3%), energy baseline calibration and others (see 
Ref.[?, ?,|5|) and statistical errors a s t a t ~ y/S. The nominal periods of data taking 
that we consider, 1 and 3 yrs, are generical representative cases where systematical 
or statistical errors are taken as predominant. 

The solar neutrino contribution can be written in the following way: 

2 2,2,2 { a\ 

Xq = Xgiob + XSK +XSNO- ( 4 ) 

glob 



The function xl\ b correspond to the total event rates measured at the Homes 



take experiment [29| and at the gallium experiments SAGE [30, 31 1, GNO |33] and 



GALLEX |33J. We follow closely the definition used in previous works (sec Rcf. ]14] 



8 



for definitions and Table (1) in Ref.[14| ' for an explicit list of results and other 
references) . 

The contribution to the x 2 from the SuperKamiokande data (x|k) nas been ob- 
tained by using double-binned data in energy and zenith angle (see table 2 in Ref. [43] 



and also Ref. ) : 8 energy bins of variable width and 7 zenith angle bins which in- 
clude the day bin and 6 night ones. The definition is given by: 

_ /^Rth _ raexpxt / 2 ,2 \~1 f n Tfth _ R cxp 



XSK 



(aR th - R^Y Knc + 4> r ) («« th - ^ CXP )- (5) 



The theoretical and experimental R quantities are this time matrices of dimension 
8x7. The factor a is a flux normalization with respect to the value measured by SNO 
NC. The covariance quantity a is a 4-rank tensor constructed in terms of statistic 
errors, energy and zenith angle bin-correlated and uncorrelated uncertainties. The 
data and errors for individual energy bins for SK spectrum has been obtained from 



Ref. @. 



The contribution of SNO to the x is given by 

2 (a - a th \ 2 2 
XSNO = I ~ ) +Xspcc-SNO ( 6 ) 

The presence of the first term is due to the introduction in our analysis of the flux 
normalization factor a with respect to the SNO NC flux, whose central and error 



values are given in Table I of Ref.[14|. The second term in formula y corresponds to: 



X 2 P ec-SNO = ^(a^-^^^aLt + ^st)" 1 ^-^ 15 ), (7) 



1.1 1 



where the day and night R vectors of dimension 17 are made up by the values of the 
total (NC+CC+ES) SNO signal for the different bins of the spectrum. The statistical 
contribution to the covariance matrix, <7 s t a ti is obtained directly from the SNO data. 
The part of the matrix related to the systematical errors has been computed by us 
studying the influence on the response function of the different sources of correlated 
and uncorrelated errors reported by SNO collaboration (see table II of Ref. ||]), we 
assume full correlation or full anticorrelation according to each source. 

To test a particular oscillation hypothesis against the parameters of the best fit 
and obtain allowed regions in parameter space we perform a minimization of the 
three dimensional function x 2 (Am 2 , tan 2 8, a). For a = a m - m , a given point in the 
oscillation parameter space is allowed if the globally subtracted quantity fulfills the 
condition Ax 2 = x 2 (Am 2 ,6 l ) - xL < xl(CL). Where x 2 = 3 (90%, 95%, ...) are the 
quantiles for three degrees of freedom. 

In Fig.|2| we graphically show the results of this analysis. They represent exclusion 
plots including KamLAND global rates, given a hypothetical experimental global 
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signal ratio: respectively strong and medium suppression S/So = 0.3,0.6 and no 
oscillation evidence S/Sq = 1.0 for one and three years of KamLAND data taking. 
As can be seen in the figures, as the KamLAND experimental signal decreases, the 
LMA region is singled out. The periodic shape in Am 2 of the 90 % C.L. (red regions) 
which becomes apparent in the three-year plot is due to the periodicity of the response 
function: in order to distinguish among these different equally-likely solutions, one 
must analyze the energy spectrum, this will done in the main analysis to be presented 
in the next sections. Obviously, only if KamLAND sees some oscillation signal (i.e. 
Si/So « 1.0 ) does the LMA region become the only solution. If we consider a 
hypothetical signal closer to 1.0 than 0.3, we see that the LOW region survives, 
although it is less favored. 



4.2 Analysis B 

Here we use the expected binned KamLAND signal for some benchmark, arbitrar- 
ily chosen, points in parameter space that we show in table ([!]). For any of these 
points we obtain the expected spectrum after 1 or 3 years of data taking under the 
"standard" conditions described above. Next we perform an standard x 2 analysis 
introducing statistical and assumed systematics errors including the evidence of the 
up-date solar experiments (CL,GA,SK and SNO). 

In the present study the total x 2 value is given by the sum of two distinct contri- 
butions, that is the one coming from all the solar neutrino data and the contribution 
of the KamLAND experiment: 

X 2 = X© + X 2 S pec,KL- ( 8 ) 

The contribution of the solar neutrino experiments x% ls described in detail in the 
previous section. The contribution of the KamLAND experiment is now as follows: 

Xspec.KL = (-R th '° - -R th ) T [ptnc + ^Iot^) ^ {R^ ~ R th ) (9) 

Note that the addition to x 2 °f a constant term Nd f has no practical importance 
for the main purpose of this kind of analysis: the determination of exclusion regions 
proceeds from the minimum subtracted quantity x 2 ~ Xmin- The are length 12 
vectors containing the binned spectrum (0.5 MeV bins ranging from 2 to 8 MeV) 
normalized to the non-oscillation expectations. Theoretical vectors are a function 
of the oscillation parameters: R th = R th (Am 2 ,9). The "experimental" vectors are 
defined in similar way, for any of the benchmark points (A?tiq,6?o) we have i? th '° = 
.R th >°(Am2,6> ). 

We generate acceptance contours (at 90,95 and 99 % CL) in the (Am 2 , tan 9 ) plane 
in a similar manner as explained in the previous section. For the sake of comparison 
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we have also obtained exclusion regions derived from the consideration of the x% e c kl 
alone. 

In Figs.(||) we graphically show the results of this analysis for a selection of points 
and for three years of data taking restricting ourselves to the LMA region of the pa- 
rameter space where, as we have noted before, the KamLAND spectrum information 
is specially significative. In each plot allowed regions corresponding to different start- 
ing points are superimposed, every region is distinguished with a label. The position 
of initial points is labeled with solid stars. 

The first case, study of the KL spectrum alone (xi pec kl) ls represented by the 
Fig. (H] left). The allowed parameter space corresponding to each particular point is 
formed by a number of, highly degenerated, disconnected regions symmetric with 
respect the line tan# = 1. These regions can extend very far from the initial point 
specially in terms of Am 2 but also in some occasions in terms of tan 2 9. For example 
the point "A" located at (Am 2 = 5.7 x 10~ 4 ,tan 2 6> = 0.38) gives rise to two sets of 
thin regions situated respectively at Am 2 ~ 10~ 3 ,10 -4 and a third region situated 
at Am 2 ~ 10~ 5 which practically covers the full range tan 2 9 ~ 0.1 — 10. A similar 
behavior is observed for point "B". Of course this situation is not very favorable 
for the future phenomenologist trying to extract conclusions from the KamLAND 
data. A much comfortable situation is found for points nearer the center of the LMA 
region. Note how the regions corresponding to the points "D,E" and specially "F" 
only extend very gently around the initial location. 

The results of the full analysis are summarized in Fig.(|3|, right). The position 
of the minima of x 2 ; marked in the plot with crosses, is practically identical to 
the position of the initial points except in some case where the difference is not 
significant anyway. The general effect of the inclusion of the solar evidence in the 
X 2 is the breaking of the symmetry in tan 2 9 as expected and the general reduction 
of the number of disconnected regions corresponding to each point. Note however 
that the point "A" still gives rise to a small allowed region situated nearly one order 
of magnitude smaller in Am 2 . The "B" region is shrinked near its initial location 
as happens to the rest of points. The conclusion to be drawed from these plots is 
that KamLAND together with the rest of solar experiments will be able to resolve 
the neutrino mixing parameters with a precision of <51ogAm 2 ~ ±0.1 practically 
everywhere. However, for values of Am 2 > 10~ 4 the problem of the coexistence of 
multiple regions with similar statistical significance will still be present. 

5 Summary and Conclusions 

We have analyzed the present experimental situation of our knowledge of the neutrino 
mixing parameters in the region of the parameter space that is relevant for solar 
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neutrinos and we have studied in detail how this knowledge should improve with the 
forthcoming reactor experiment KamLAND. 

In this work we have presented in some detail the characteristic of KamLAND ex- 
periment including antineutrino reactor fluxes and absortion cross sections for which 
we have given updated values. We find that present theoretical errors in this cross 
secion are negligible. 

We have studied the expected KamLAND spectrum for different possible ( "bench- 
mark points" ) values of the mixing parameters selected inside the LMA region. The 
shape of the spectrum shows a significant dependence on the values of the mixing 
parameters. The spectrum distortion caused by oscillation have been characterized 
by the first moment of the positron energy distribution. The results confirm that 
KamLAND is very sensitive to neutrino oscillation in the LMA region. In particular 
the value of the first moment changes very much for values of Am 2 varying inside 
the region Am 2 ~ 10~ 5 — 10~ eV 2 . The dependence on the value of the mixing angle 
is also evident, even if somehow milder. In the LOW and SMA region, instead, the 
value of the moment is essentially constant. We have also verified that the result 
doesn't depend in a significant way on the choice of the bin size. 

In order to investigate the discrimination power of KamLAND, we have selected 
different points ("benchmark points") in the LMA region and studied which infor- 
mation KamLAND will be able to give after 1-3 years of running. We have included 
a full modelling of statistical and systematic uncertainties. The regions selected by 
KamLAND alone, all symmetric with respect to tan 2 9 = 1, have a large spreadth in 
the mixing angle. The experiment should have, instead, a much higher sensitivity to 
the mass difference parameter, especially for values of Am 2 lying in the central part 
of the LMA region. The situation would be more complicate for values closer to the 
border of the LMA region or beyond already in the HLMA region (i.e. Am 2 < 2xl0 5 
and Am 2 > 8-9 x 1CT 5 or tan 2 far from 0.5). In this case KamLAND, with or 
without solar evidence, will be able only to select multiple regions in the parameter 
space, in the sense that different possible values of the parameters would produce the 
same signal. 

We have performed a similar analysis adding to the information from KamLAND 
would-be signal all the evidence already exististing on solar neutrinos. By using a x 2 
analysis, we have produced exclusion plots. KamLAND will help to select the values 
of the mixing parameters especially in the case in which the solutions lies in the LMA 
region. If, instead, one moves towards values of the KamLAND signal closer to the 
no oscillation value (that at present seems to be strongly disfavored by the other 
experiments) the absolute x 2 minimum moves from the LMA to the LOW solution 
and one is left with small allowed regions not only in the LOW, but also in the SMA 
region. In summary, KamLAND together with the rest of solar experiments will 
be able to resolve the neutrino mixing parameters with a high precision practically 
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everywhere. However, for values of Am 2 > 10~ 4 (graphically emphatized by the two 
benchmark points labeled "A" and "B" in the plots) the problem of the coexistence 
of multiple regions with similar statistical significance will still be present. 
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Label 


tan 2 6> 


Am 2 (eV 2 ) 


A 


0.38 


5.70 x 10" 4 


B 


0.60 


2.04 x 10" 4 


C 


0.50 


1.01 x 10" 4 


D 


0.50 


3.60 x 10- 5 


E 


0.56 


2.37 x 10^ 5 


F 


0.99 


5.07 x 10~ 5 



Table 1: Benchmark points in the 2-neutrino parameter plane used in the analysis. 



3.8 




Figure 1: (Left) Different kamLAND spectra for points in the LMA region. (Right) Kam- 
LAND moments distribution as a function of the oscillation parameters. The No-oscillation 
mean-spectrum energy is given by the horizontal black line. 
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Figure 2: Exclusion plots including KamLAND global rates (Analysis A), given a hypothetical 
experimental global signal ratio: respectively S/Sq = 0.3,0.6, and no oscillation evidence S/Sq = 
1.0. Statistical and assumed systematics (~ 5%) errors are included, left (right) panel refers 
to one (three) year of KamLAND data taking. The colored areas are allowed at 90, 95, 99 and 
99.7% CL relative to the absolute minimum. The region above the upper thick line is excluded 
by the reactor experiments Q. 
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Figure 3: Allowed areas in the two neutrino parameter space after 3 years of data taking in 
KamLAND (Analysis B). Allowed regions belonging to the same point are labeled with the 
corresponding letter, the position of the point itself is labeled with a solid star (see table [j]). 
The colored lines separate allowed regions at 90, 95, 99 and 99.7% CL relative to the absolute 
minimum. (Left) Results with the KamLAND spectrum alone. (Right) KamLAND spectrum 
plus solar (CL,GA,SK,SNO) evidence. Crosses are situated in the position of the \ 2 minima. 
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